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ABSTRACT 

Nearby starburst galaxies liave consistently shown anomalous Lya/H/3 ratios. By re-analysing the published 
lUE/optical observations, we show that most starbursts present a normal Lya emission, consistent with case 
B recombination theory, provided extinction laws appropriate to their metallicities are used. This implies that 
extinction is more important than multiple resonant scattering effects. The anomalous emission and absorption 
lines present in a few remaining galaxies are simply explained if they are observed in the post-burst phase, 
between about 10^ and 10® yrs after the start of the burst. We use updated stellar population synthesis models 
to show that anomalous ratios are produced by the aging of stellar populations, since the underlying stellar Lya 
line is important in the cooler massive stars. The inferred low-duty cycle of massive star formation accounts 
naturally for the failure to detect large numbers of Lya-emitting galaxies in deep surveys and at high redshift. 
Some testable predictions of the proposed scenario are also discussed. 
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1 INTRODUCTION 



In any galaxy formation scenario young galaxies are 
predicted to present a strong Lya emission line, essentially 
associated with the cooling of the gas and the subsequent 
formation of stars. The details of the emission depend on 
the assumed cooling and star formation rates during the 
collapse and infallof gas (Partridge & Peebles 1967, Kauf- 
man 1976, Meier 1976, Rees 1988). Recent simulations of 
an inhomogeneous dissipational collapse (Baron & White 
1987) also indicate a strong Lya line. 

Yet despite the large number of observations aimed at 
detecting such an emission, no strong Lya sources have 
been found at high redshift that can be classified as pro- 
togalaxies (Davis 1980, Koo 1986, Pritchet & Hartwick 
1987, Lowenthal et al. 1990, De Propis et al. 1993); both 
distant quasars and radio-galaxies have a Lya emission 
probably dominated by non-thermal processes. A possi- 
ble exception are the recently observed Lya emission lines 
associated with damped Lya systems (Hunstead, Pettini 
k Fletcher 1990, Lowenthal et al. 1991, Wolfe et al. 1992, 
M0ller k Warren 1993). 

One way to understand the conditions of Lya emis- 
sion at high redshift is to study the properties of nearby 
starburst galaxies. Several observations have consistently 
shown that their Lya emission is quite peculiar, being very 
weak, absent or even in absorption (Meier k Terlevich 
1981 [MT], Hartmann et al. 1984 [H84], Deharveng, Jou- 
bert k Kunth 1985 [DJK], Hartmann et al. 1988 [H88]). 
It was thus argued that their corresponding Lya/II/3 ra- 
tios are inconsistent with case B recombination, rendering 
difficult the interpretation of their spectra. The current 
explanation is that the multiple resonant scattering with 
Hi atoms increases the path length of Lya photons, and 
accordingly the probability of absorption by dust (Hum- 
mer k Kunasz 1980). For this to hold, however, the opti- 
cal depth of dust at Lya must be relatively large and an 
homogeneous Hi envelope with unit covering factor must 
surround the sources. In this paper we show that this is 
not the case, and we suggest an alternative mechanism. 

2 THE LYa/H/3 RATIO IN STARBURSTS 

We consider the entire sample of nearby starburst 
galaxies with lUE and optical observations (MT,H84,DJK, 
H88, Margon et al. 1988, Calzetti k Kinney 1992, Ter- 
levich et al. 1993 [T93]). In order to disentangle extinc- 
tion corrections from other effects, we study the depen- 
dence of the extinction-corrected \jy a /Mj3 emissivity ratio 
on the observed (uncorrected) Ha/H/3 fiux ratio (Figure 
1). The horizontal axis in this diagnostic diagram is thus 
a function of the optical depth of dust, while the ver- 
tical axis measures the deviation from standard recom- 
bination theory. If resonant scattering effects are unim- 
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Figure 1: Ratio of de-reddened Lya/H/3 emissivities as a 
function of the observed Ha/H/3 fiux ratio. The extinc- 
tion law appropriate to the metallicity of each galaxy is 
indicated. The filled triangle and circle represent the po- 
sitions of Mrk 357 and Mrk 496, respectively, at small 
apertures. Measured with a large (lUE) aperture, their 
positions (open symbols) agree with standard recombina- 
tion theory. 

portant, galaxies should lie between the dotted lines in 
Figure 1. To derive the corrected Lya/H/3 ratio the ex- 
tinction law appropriate to the metallicity of the galaxy 
must be used. We adopt characteristic mean extinction 
laws corresponding to Galactic, LMC and SMC metallic- 
ities (Seaton 1979, Fitzpatrick 1986, Bouchet et al. 1985) 
to correct the ratios of each galaxy. 

As Figure 1 shows, most of the galaxies do indeed lie 
close to the range predicted by standard recombination 
theory, even when the dust optical depth is large (large 
Balmer decrement) and regardless of the Hi column den- 
sity present in each galaxy. Thus most galaxies have 
a normal Lya emission and their Lya/H/3 ratios agree 
with case B recombination, provided appropriate extinc- 
tion laws are used. This in turn implies that resonant 
scattering effects are negligible, a fact that can easily be 
interpreted within the framework of inhomogeneous, mul- 
tiphase models for the interstellar medium of starbursts 
(Neufeld 1991, Valls-Gabaud 1993 [VG93]) in which Lya 
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photons escape through the hot, ionised phase. The an- 
ticorrelation between the Lya/H/3 ratio and the oxygen 
abundance that was found previously (DJK,MT,H88,T93) 
is not confirmed by our analysis, since galaxies in the en- 
tire range of metallicities, from Galactic (0833+652, Mrk 
496) to SMC (Pox 120, Tol 41), present normal Lya/H/3 
ratios. 

Some galaxies in the upper part of Figure 1 might have 
been overcorrected, in particular Tol 41, but this may 
be due to several effects. For instance the Balmer decre- 
ment can overestimate the extinction (as measured by the 
thermal radio continuum, e.g. Lequeux et al. 1981). Fur- 
thermore the extinction in the UV might be lower than 
in the visible (Fanelli et al. 1988). We expect that the 
high energy density contained in the radiation field of the 
starburst environment efficiently destroys the PAHs, as 
observed for example in the LMC (Sauvage, Thuan & Vi- 
groux 1990). The extinction at Lya must then decrease 
since PAHs make up an important part of the UV extinc- 
tion (Desert, Boulanger & Puget 1990). Also, since the 
oxygen abundance is enhanced during a starburst an over- 
estimation of the actual metallicity is produced (Gilmore 
& Wyse 1991) and may change the correction too. 

A third of the sample lies in the lower part of Figure 
1, and these galaxies (Pox 124, Mrk 26, Tol 1247-232 and 
C0840-I-120) are truly anomalous, since both their dust 
optical depth and their Lya emission are small. Several 
effects can explain the discrepancy. The most important 
one is possibly due to the different apertures used for the 
lUE spectra and the optical ones. Calzetti & Kinney 
(1992) have reobserved in the optical three galaxies us- 
ing the same (large) aperture as lUE, and found that the 
Balmer decrements were significantly larger than those 
measured using standard small apertures. In Figure 1 the 
positions of Mrk 357 (filled triangle) and Mrk 496 (filled 
circle), as derived from spectra taken with small aper- 
tures are indicated. Small apertures seem to underesti- 
mate the integrated Balmer decrement, and accordingly 
the Lya/H/3 ratio, so it seems that aperture effects may 
account for the apparent anomalous ratios. Another pos- 
sibility is of course the destruction by enhanced scattering 
in a uniform Hi envelope and absorption by dust. Yet why 
this process should be efficient in anomalous galaxies but 
not in normal ones is unclear. Moreover, this effect would 
require an improbably large Hi column density of 10^^ 
cm~^ for Pox 124. 

Even if these effects may account for the remaining 
anomalous galaxies, neither of them explains the presence 
of Lya lines in absorption detected in Mrk 309, Mrk 347 
and Mrk 499 (H88). Another mechanism must necessarily 
play an important role. We suggest below that this is the 
aging of the stellar populations. 



3 LYa EQUIVALENT WIDTHS AND AGING 

In principle, if an homogeneous Hi layer covers most 
of the Hll regions, the equivalent width (EW) of the Lya 
line in a galaxy spectrum is difficult to interpret, since line 
photons will be affected by multiple resonant scatterings, 
while continuum photons will not. However, this effect 
is negligible in nearby starbursts (§2), and thus the Lya 
EW can be used to infer the properties of the underlying 
stellar populations. 

An important and previously unaccounted-for effect 
which strongly reduces the Lya emission is the fact that 
late B stars have very large Lya equivalent widths in ab- 
sorption, as Figure 2 shows (Valls-Gabaud 1991). The 
earlier B and O stars (Teff > 22,000 K) have small EWs, 
but the later ones show an exponential increase with de- 
creasing temperature. This trend is clearly confirmed by 
OAO-2 (Savage & Panek 1974) and Copernicus observa- 
tions (Vader et al. 1977), in addition to being consistent 
with predictions of NLTE models (Fig. 2). In the 30,000 
K < Teff < 10, 000 K range, the trend for main-sequence 
stars is 
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where absorption EW is denoted with a minus sign. 
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Figure 2: Equivalent width of the stellar Lya absorption 
line as a function of effective temperature. The continuous 
line gives the theoretical predictions from NLTE models, 
while the dotted line is an empirical trend (Eq. 1). 
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No observations were found for early A type stars, ex- 
cept for Vega which has about 100 A (Praderie 1981) and 
confirms the trend given by Eq. (1). Although lower 
temperature (F, G and M) stars may have a Lya emis- 
sion produced by the chromosphere (see e.g. Landsman 
& Simon 1993), their contribution to the integrated UV 
fiux is negligible for ages less than 10® yrs. Massive OB 
supergiants might develop a P Cygni profile in the Lya 
line with a net EW in emission of less than 1 A (Mihalas 
& Hummer 1974) although the details depend strongly on 
the wind properties (Sellmaier et al. 1993). In any case, 
this shorter lived phase does not affect the predictions for 
integrated spectra for ages less than 10® yrs which will be 
dominated (at Lya) by main sequence stars. 
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Figure 3: Evolution of the Lya equivalent width for a 
burst with the three indicated IMFs and a constant star 
formation rate (miom = O.IMq in all cases), (a) Total 
equivalent width (emission-|- absorption). After ~6 10® 
yrs, the EW is independent of the IMF. (b) Integrated 
absorption equivalent width for the same cases. 



We have included the effect of the stellar Lya line in 
the updated stellar population synthesis code of Bruzual 
and Chariot (1993) to estimate the net EW of the ac- 
tual line. Figure 3a shows the evolution of the total EW 
(Wtot(Lya)) as a function of time for a burst and a con- 
stant star formation rate. The results are very robust and 
do not depend strongly on the slope a of the initial mass 
function (IMF), nor on its upper cutoff m^p, except in the 
very early phases (t < 10® yrs). During the first 10®'® yrs, 
the EW is about 250 A, decreasing sharply afterwards so 
that by 10^ yrs the Lya line becomes negligible or indeed 
is in absorption. The time scale corresponds to the life- 
time of the later B stars which are the only significant 
contributors to the far UV at that time. Note that the 
strong increase after 10® yrs is produced by the emission 
from planetary nebulae, which is irrelevant for the galax- 
ies studied here, but may be important in other galaxies 
(Hansen et al. 1991). 

Figure 3b shows the large contribution of the underlying 
stellar population to the absorption EW (Wabs), typically 
50 A at 10® yrs in the burst case. In the constant star for- 
mation rate case it is much less important since massive 
OB stars are continuously being formed and dominate the 
UV spectrum, and their integrated Wabs(~^10 A) is al- 
ways much smaller than the emission one. The use of low 
metallicity stellar evolutionary tracks and spectra (Olof- 
sson 1989, Mas-Hesse & Kunth 1991) does not change 
these values significantly. 

In order to test whether these results are actually ob- 
served, we plot in Figure 4 the net Lya EW versus the 
(U-B) colour for the sample of nearby starbursts (crossed 
circles). Galaxies with no published (U-B) colours were 
placed at the mean and median colour of the sample, (U- 
B)m = -0.4 (open circles) and are likely to be in the range 
-0.3 to -0.6. While the constant SFR is ruled out, given 
the distribution of galaxies in this diagram, the burst case 
provides an excellent fit and indicates several interesting 
properties: 

i) in the same range of very blue colours there is a con- 
tinuous distribution of galaxies, ranging from strong Lya 
emitters to galaxies showing Lya in absorption. 

ii) the main parameter of this distribution is the age, since 
for populations younger than about 10^ yrs the Lya emis- 
sion is normal and the EW is large, while for older bursts 
the Lya line is smaller or becomes in absorption, indepen- 
dently of the details of the IMF. 

The aging of the stellar populations in a burst is thus a 
simple unifying mechanism that naturally accounts for the 
presence of normal Lya emitters, galaxies with 'anoma- 
lous' Lya/H/3 ratios, and Lya lines observed in absorp- 
tion. If a galaxy is observed in the 'active' burst phase, 
normal Lya/H/3 ratios are observed. This is the case for 
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Figure 4: Lya equivalent width as a function of the U — B 
colour, for the same star formation rates and IMFs indi- 
cated in Fig. 3. Observations (crossed circles) and esti- 
mated uncertainties are marked. Open circles represent 
galaxies with measured Lya EWs, but with no published 
colours (see text). The predicted position of flat spectrum 
galaxies is also indicated. 



most of the galaxies in the sample. Since the lines are un- 
resolved, the measured flux scales as (1 — | VFabs l/M^tot) and 
decreases with time (Fig. 3) so that seemingly 'anoma- 
lous' Lya/H/9 ratios appear if the galaxy is observed in 
the 'post-burst' phase, about 10^ yrs after the start of the 
burst. When the burst is observed at an older age, be- 
tween a few 10^ and 10® yrs, the Lya line is in absorption 
(Fig. 4). 

The Lya line is far more sensitive to the aging of the 
stellar populations than the Balmer or the [Oll]A3727, 
[Olll]AA4959,5007 lines since the stellar absorption line 
is much larger for Lya than for the latter lines. There- 
fore, unlike the star formation rate obtained through the 
Ha emission (Kennicutt 1983), for which the stellar ab- 
sorption line is always relatively small (e.g. Diaz 1988), 
a star formation rate derived from the Lya line is highly 
unreliable, essentially because SFR(Lya) cx EW(Lya) x 



SFR(UVcont). A rate of star formation derived from the 
UV continuum flux or Balmer emission lines is more suit- 
able since it is less biased. A contribution by supernovae 
remnants or shocks is not required to explain the Lya line 
in these starbursts. 

4 IMPLICATIONS FOR HIGH REDSHIFT GALAXIES 

Ultra-deep imaging surveys (eg. Tyson 1988, Cowie 
et al. 1988) have revealed a population of possible 'flat- 
spectrum' galaxies. Our proposed scenario predicts that 
for those galaxies with flat rest-frame spectra (U-B~- 
0.9), a conspicuous Lya emission should be present with 
large EWs of the order of 250 A (see Figure 4). Yet we 
expect that only a small fraction will be visible due to 
the short lived nature of the flat spectrum phase in U-B. 
Furthermore if high redshift radio galaxies have similar 
spectra (eg. Lilly 1989), the maximum EW expected from 
the associated young stellar population cannot exceed 300 
A, which is close to the observed values, although a non- 
thermal contribution by the nucleus or by shocks cannot 
be excluded. 

In addition, this scenario easily explains the lack of 
success of deep imaging surveys in detecting primordial 
galaxies through their Lya emission. Since any SFR can 
be decomposed into a series of short bursts (e.g. Searle, 
Sargent & Bagnuolo 1973), the Lya emission from a 
galaxy is a strongly varying function of time. If Ni, is 
the number of bursts per galaxy in a Hubble time, and 
ti, the duration of the burst in Lya, the fraction /j of 
galaxies that are observable in Lya emission is 



fb = Niti/tH ~ 10-3 h Ni 



(2) 



The number of bursts Ni, must satisfy the constraint that 
the total mass in stars formed in a Hubble time is the 
same as that produced by the average SFR in that period, 
hence 



50 h- 



SFRt/ < SFR > 
20 



lO^yrs 



(3) 



Then Eq. (2) implies that in the deep volumes sampled 
by the present surveys, only about 5% of the galaxies 
will be visible in Lya. The constraints set by Lya sur- 
veys (Pritchet & Hartwick 1990, Lowenthal et al. 1990, 
De Propis et al. 1993) on the theoretical scenarios of 
galaxy evolution are therefore much less stringent, by at 
least one order of magnitude, on the number density of 
Lya emitters at a given flux. We also note that this 
argument explains the low rate of Lya emission detec- 
tions in damped Lya systems (VG93). Ongoing imaging 
(Thompson, Djorgovski & Trauger 1992) and spectral sur- 
veys (Cowie et al. 1988) can provide a further test of this 
mechanism by setting limits on the duty cycle. 
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5 DISCUSSION AND CONCLUSIONS 
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